Three-dimensional shear in granular flow 
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The evolution of granular shear flow is investigated as a function of height in a split-bottom 
Couette cell. Using particle tracking, magnetic-resonance imaging, and large-scale simulations we 
find a transition in the nature of the shear as a characteristic height H* is exceeded. Below H* 
there is a central stationary core; above H* we observe the onset of additional axial shear associated 
with torsional failure. Radial and axial shear profiles are qualitatively different: the radial extent is 
wide and increases with height while the axial width remains narrow and fixed. 

PACS numbers: 45.70.Mg, 45.70.-n, 83.50.Ax 



Shear bands in dense granular materials are localized 
regions of large velocity gradients; they are the antithesis 
of the broad uniform flows seen in slowly-sheared New- 
tonian fluids P, 0, S 0) 0) 0] ■ Until recently it was gen- 
erally assumed that all granular shear bands were nar- 
row. However, in 2003 Fenistein et aL0] discovered that 
in modified Couette cells granular shear bands can be 
made arbitrarily broad. In this geometry, the bottom 
of a cylindrical container is split at radius r = Rg and 
shear is produced by rotating both the outer ring and the 
cylindrical boundary of the container while keeping the 
central disk (r < Rs) stationary. For very shallow packs, 
the shear band measured at the top surface is narrow and 
located a.t r — Rg so that the inner region directly above 
the central disk is stationary while the remaining part ro- 
tates as a solid. As the filling height of the material, H, 
increases, the shear band increases in radial width and 
moves toward the cylinder axis. For sufficiently large H, 
the shear band overlaps the axis at r = and one miglit 
expect qualitatively new behavior. Indeed, Unger et al.^ 
predicted that the shape of the boundary between mov- 
ing and stationary material would undergo a first-order 
transition as H is increased past a threshold value H*: 
the shearing region which for H < H* is open at the 
top and intersects the free surface abruptly collapses to 
a closed cupola completely buried inside the bulk. 

Previous experiments focused primarily on the surface 
flows in shallow containers and left unexplored many 
questions about the shape and evolution of the shear 
profiles for large H. Here, we combine magnetic reso- 
nance imaging (MRI) and high-speed video observations 
with large-scale simulations to explore shear flow both 
for shallow and tall packs. In addition to monitoring the 
evolution of the flow profiles in the radial direction, we 
also examine shear in the vertical direction. Instead of 
a first order collapse of the shear zone as proposed by 
Unger et a/.Q, we find that above H* ~ Q.GRg, the inner 
core of immobile material disappears gradually as shear 
along the central axis of the cylinder sets in. 
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FIG. I: Surface fiow as function of filling height, H. (a) 
Angular velocity, u>o at the cell axis, (b) Center of the shear 
band, Rc. Solid line is a gaussian fit to guide the eye; dotted 
line is the fitting function of Fenistein et al. 0; dashed line 
is the theoretical result 0. Vertical dash-dot line indicates 
H = O.GRs- Inset: Schematic of split-bottom Couette cell. 



Our setup is similar to that of Fenistein et aL0] ex- 
cept that we rotate the inner disk instead of the outer 
ring and cylinder (Fig^D inset). In the absence of iner- 
tial effects, this makes no difference to the results. For 
surface observations with high-speed video we use a cylin- 
drical cell with radius Rout = 72.25'mm whose bottom 
is split at Rs = 55.5mm. For the MRI experiments a 
cell with Rout = 43.9mTO and Rs = 32.5mm is used. 
The inner disk is rotated at angular velocities, fl, in the 
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shear-rate-independent regime. For the data presented 
here = 1.96rad/s and fl = l.Olrad/s for the surface 
and MRI experiments, respectively. We fill the cell with 
granular material to a total height h = H, measured from 
the cell bottom. A layer of grains glued to the cell walls 
and bottom assures controlled friction at the boundaries. 
For the surface flow measurements we used spherical mus- 
tard seeds {d = 1.9mm) and tracked their motion with 
high-speed video at frame rates ranging from 250s^^ to 
0.027s~^. To follow the motion of particles inside the 
pack using MRI, we use a mixture of poppy and rajagara 
seeds. Rajagara seeds are more spherical than poppy 
seeds(Fig|2l3 inset), but have nearly the same average di- 
ameter ((i=0.85mm) and the same density (p=l.lg/cm'^). 
Poppy seeds contain more oil than rajagara seeds pro- 
viding a clear contrast in MRI signal which allows for 
particle tracking (Fig|2K inset). 

The simulations are carried out with a discrete ele- 
ment method in which grains interact only upon con- 
tact through assumed point forces in normal and tan- 
gential directions and elastic tangential displacements 
are truncated as necessary to satisfy the Coulomb cri- 
teria at the contact. Details of the specific implementa- 
tion can be found elsewhere 0. We use mono-disperse 
hertzian spheres with a layer of frozen particles at the 
bottom. The relevant parameters describing the ma- 
terial properties of the spheres are the normal stiffness 
kn — 2 X lO^mg/d, the tangential stiffness kt — 2/7fc„, 
the normal and the tangential viscous damping coeffi- 
cients 7„ — 50-\/ g/d, 7t — 0, and the particle and wall 
coefficients of friction /i = 0.5, where d and m are the di- 
ameter and mass of spheres and g is the gravity acceler- 
ation. We have checked that changing the coefficients of 
friction does not qualitatively change the observed behav- 
ior. The cell dimensions used for the simulation match 
those used for the surface measurement experiments. We 
choose n = 0.014y^ ~ 1.39md/s. 

Evidence for a change in flow behavior with increasing 
filling height, H, can already be found by tracking par- 
ticle motion at the free top surface (Fig^. The angular 
velocity at the center of the cell, cuq — uj{r — O^h = H) 
is independent of H for shallow packs but begins to de- 
crease beyond H* ~ 0.6i?s. At a similar height Fenistein 
et al. found that the radial shear profile at the surface 
begins to deviate from the universal error-function shape 
describing the profile for shallow packs 0]. Dramatic 
deviations are also seen in the evolution of the center 
position of the shear band, Rc{H), with filling height 
(Fig^). Our results for Rc, defined as the radial posi- 
tion where the shear rate has its maximum, are consistent 
with previous experiments as well as with models jsj 
in the shallow pack regime. However, while tuo decreases 
for H > H*, the shear zone does not disappear at the 
surface as predicted by the theory (dashed line). Instead, 
in both experiments and simulations Rc asymptotically 
approaches a nonzero value. 
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FIG. 2: Angular velocity profiles uj{r) at different vertical 
positions h for (a) H = OARs and (b) H = O.SSRs- Symbols 
for MRI experiments: (a) h = 0.21// (D), h = 0.43// (o) 
and h = 0.65// (V); (b) h = 0.10// (□), h = 0.21// (o), 
h = 0.31// (A), 0.42// (V) and h = 0.73// (o). Lines for 
simulations; (a) from right to left h = 0.19//, h = 0.50// and 
h = 0.85//; (b) from top to bottom h = 0.11//, h = 0.23//, 
h = 0.30//, h = 0.43// and h = 0.70//. (a) Inset: MRI image 
for one layer inside the bulk. Bright spots are poppy seeds, 
dark background are rajagara seeds; poppy seeds glued to the 
wall of the cell show up as a bright circle, (b) Inset: Optical 
micrographs of rajagara (top) and poppy seeds (bottom). 



These results imply that, beyond H*, velocity gradi- 
ents must also exist in the vertical, axial direction near 
the center of the cell. Our MRI experiments and simu- 
lations explore this shear flow inside the bulk. We pre- 
pare our MRI samples by mixing 5% (by volume) poppy 
seeds (MRI positive seeds) uniformly with rajagara seeds 
(Fig|21i inset). Images before rotation and after an in- 
terval of rotation are taken. By performing a cross- 
correlation of the two images as a function of radius, we 
obtain velocity profiles a;(r) as a function of h (FigjSJ. 
This method enables us to measure velocities with sev- 
eral orders of magnitude difference (for details, see ^H). 
For both experiments and simulations care was taken to 
assure that the systems are in the steady state by rotat- 
ing long enough before any measurement were performed. 
As an additional check we made sure that stopping and 
restarting the system did not change the velocity profile, 
which is consistent with the previous study [lOl |. 

For H < H* , MRI and simulations show an inner core 
at the center of the pack which rotates as a solid along 
with the rotating bottom. We plot the profiles uj{r) in 
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FIG. 3: (color). Contours of constant angular velocity, for different filling height H. Upper panel: MRI experiments. 

Lo/Q = 0.84(0), 0.24(B), 2.4 x 10"^(o), 2.4 x fO"^(A), 2.4 x 10"''(V). Dashed lines indicate H and dotted lines are guides to 
the eye. Lower panel: simulations. Color is used to identify velocity ranges. Dark red: ui/Q G [0.84, 1]; orange: [0.24,0.84]; 
yellow: [2.4 X 10"^0.24]; green: [2.4 x f 0"^ 2.4 x 10"^]; blue: [2.4 x 10"*, 2.4 x 10"^]; white: [0,2.4 x 10"'']. 



FigEK- There is no slip between different layers near the 
center of the cell, and the profiles are fit well by an error 
function. However, when wc increase H above H*, axial 
slip occurs in both MRI experiments and simulations: 
the bottom layer rotates at the same rate as the bottom 
disk, while the layer near the surface hardly rotates at all 
(Figl^b) 'ill- Thus the decrease in surface flow velocity 
for H > H* is caused by shear between horizontal layers 
inside the bulk. For H < H*, an inner core exists at the 
center of the cell which rotates with the inner disk with 
^^\r=o = 0; while for H > H\ ^^\r=, ^ 0. 

To visualize the resulting shear profiles, we plot cross- 
sections of the system with contours of constant angular 
velocity (Fig|2|)- For H well above H* , the high- velocity 
contours close into dome shapes (FigOli), which gradu- 
ally open as H decreases. The contours for smaller ve- 
locities open up earlier and eventually touch the surface 
(Fig|3l3, c). When H < H*, all velocity contours touch 
the surface as a solid inner core forms for r < Rg (Fig|3^). 
The lower panel shows the corresponding simulation re- 
sults in a color gradient. One can see how the motion at 
different h correlates with the motion of the bottom disk. 

Having access to the full shear profiles allows us to 
address the question whether the same length scale de- 
scribes the shear in axial and radial directions. Fig0| 
inset a, b show uj{h) measured along the axis of the cell 
(r = 0) for different heights H. We note that the axial 
velocity decays in an exponential fashion before the data 
plateau at a level that extends to the surface. All MRI 
and simulation data can be fit consistently to a gaus- 
sian form, u){h)/Vl = a + (1 — a)cxp[— /i^/(2(7^)], where 
a is an i?-dependent offset indicating the angular ve- 
locity at the top surface and a is the axial shear-band 
width. The main panel in Fig 0] demonstrates the consis- 
tency of fitting a gaussian to the axial shear profile, col- 
lapsing data from MRI measurements and simulations. 
Near H* , the offset, a, approaches unity in an exponen- 
tial manner (Fig[Sl inset). Although all data show the 
same exponential behavior, the offsets from simulations 
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FIG. 4: Scaled angular velocity {(lo/Q) — a)/(l — a) for dif- 
ferent heights h/a. Insets show the corresponding unsealed 
data from (a) MRI and (b) simulations. The main panel 
uses the same symbols as the insets. Solid line is a gaussian 
exp(-/i^/(2cr^)) with a = O.lSi?^. Inset (a): H = 0.97iis(B), 
H = 0.887?, (.), H = 0.78ii,(A), H = 0.74i?,(T) and 
H = 0.63i?,(*). Inset (b): H = 1.02i?,(n), H = 0.90i?,(o), 
H = 0.78i?.(A), H = 0.72i?,(v) and H = 0.66i?.(0). Solid 
lines are gaussian fits introduced in the text. 



and surface measurement with mono-disperse spherical 
grains are larger than those from the MRI experiments 
with poppy and rajagara seeds jTJj- Extrapolating each 
set of data to cj/O = 1 suggests that the onset of axial 
shear begins at H* = 0.60 ± 0.02i?,. 

The data in FigEl together with those from Fig^, 
demonstrate that the transition in shearing behavior at 
H* occurs in a continuous manner. This differs from 
the model by Unger et al. ^ which predicts a first-order 
transition at H* = 0.7Rs. While the Unger model, based 
on the idea of minimum dissipation of energy, includes 
the essential elements for the transition in the shape of 
the shear band, a key aspect not considered is the axial 
slip below H*. The necessity for such slip emerges from 
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FIG. 5: Width of the shear profiles, ct, for different filling 
heights H. Data show a both along the axial direction in the 
bulk (■ and □ from MRI experiments using different methods 
[l^ . • from simulations) and along the radial direction in the 
surface layer (x). Dashed line is cr = O.fS-Rs. Inset: offset a 
vs. H from MRI (■, □) and simulations (•). Surface data 
from FigQare shown for comparison (A). Solid line indicates 
exponential behavior. 



considering the torque balance in shallow packs. When 
the torque on the surface of the inner core exceeds the 
frictional strength at the bottom, slip will occur. The 
total torque on the surface of the inner core and the tor- 
sional strength of the contact between the bottom of the 
cell and the granular material above it can be calculated 
if the shape of the inner core hcoreix) is known. Using 
the approximate hcoreif) given in [g, it can be shown 
that this already occurs when H > 0.5Rs- Thus, already 
within the Unger model a torsional failure mode near the 
bottom should preempt any first-order transition within 
the bulk. Our data demonstrate that this torsional fail- 
ure is associated with a well-defined axial shear band at 
r = that exhibits a gaussian profile. 

Gaussian shear profiles have been observed in previous 
studies using traditional Couette cells 0- However, in 
the present geometry we find such a profile along the axial 
direction and an approximate error-function profile along 
the radial direction. The width of this axial shear zone 
appears to be independent of H, with an average value 
a = O.lSRs, implying a/d = 6.9 for the experiment and 
5.4 for the simulation (Fig|31). This is in contrast to the 
width of the radial shear profile (FiglSJ, which strongly 
depends on H and approaches as iJ ^ 0- 

The ubiquitous presence of shear bands is one of the 
crucial differences between granular materials and ordi- 
nary fluids. Understanding what gives rise to the shear 
profiles is one of the outstanding puzzles in granular dy- 
namics. The modified Couette geometry produces two 
distinct forms of shear: a radial component whose width 
grows with height as shown by Fenistein et al. , and an 
axial component with a small constant width that only 
appears when the filling height exceeds a threshold. The 



different character of the shear bands in the radial and 
axial directions shows that boundary conditions are es- 
sential for determining shear localization. Because the 
onset of axial shear in this geometry is continuous and 
controlled simply by the height of the pack, these studies 
have allowed for detailed observation of how shear can be 
initiated in the bulk. Very recent surface measurement 
of the central core procession also corroborates our 
three dimensional results. 
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